
clear all

use "data.dta"

* Descriptive statistics

centile age_cat

tab female

tab higher_ed

* Cronbach's alpha

alpha conf1 conf2 conf3 conf4

* ANOVAs

bysort contnt: sum agreem
anova agreem contnt
estat esize

bysort contnt: sum p_easte
anova p_easte contnt
estat esize

anova agreem source

anova p_easte source

anova agreem contnt##source

anova p_easte contnt##source

anova p_nwest contnt

* Difference in effects of conformity on agreement by message content

sem ///
(Conform -> conf1 conf2 conf3 conf4) ///
(p_easte <- agreem Conform) ///
(agreem <- Conform), ///
group(contnt) standardized

disp (.3503673-.3093814)
disp (.3503673-.3093814)/sqrt(.0695642^2+.0680372^2)
disp (2*(1-normal(.42121117)))

* Difference in effects of conformity on agreement by message source

sem ///
(Conform -> conf1 conf2 conf3 conf4) ///
(agreem <- Conform), ///
group(source) standardized

disp (.3391649-.3095974)
disp (.3391649-.3095974)/sqrt(.0524601^2+.0498204^2)
disp (2*(1-normal(.40868827)))

* Difference in direct effects of agreement on prejudice by message content

disp (.1513998+.2243219)
disp (.1513998+.2243219)/sqrt(.0671531^2+.0677599^2)
disp (2*(1-normal(3.9384296)))

* Note 5

sem ///
(Conform -> conf1 conf2 conf3 conf4) ///
(p_easte <- agreem Conform age_cat female higher_ed) ///
(agreem <- Conform age_cat female higher_ed), ///
group(contnt) standardized
estat teffects

disp (.094286+.14499)/sqrt(.0560577^2+.0562093^2)
disp (2*(1-normal(3.0141275)))

disp (.9073289-.4134749)/sqrt(.1388377^2+.1557535^2)
disp (2*(1-normal(2.3668951)))
